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We investigate which Jordan frame F{R) gravity can describe a Type IV singular bouncing 
cosmological evolution, with special emphasis given near the point at which the Type IV singularity 
occurs. The cosmological bounce is chosen in such a way that the bouncing point coincides exactly 
with the Type IV singularity point. The stability of the resulting F{R) gravity is examined and in 
addition, we study the Einstein frame scalar-tensor theory counterpart of the resulting Jordan frame 
F{R) gravity. Also, by assuming that the Jordan frame metric is chosen in such a way so that, when 
conformally transformed in the Einstein frame, it yields a quasi de Sitter or de Sitter Friedmann- 
Robertson-Walker metric, we study the observational indexes which turn out to be consistent with 
Planck 2015 data in the case of the Einstein frame scalar theory. Finally, we study the behavior of 
the effective equation of state corresponding to the Type IV singular bounce and after we compare 
the resulting picture with other bouncing cosmologies, we critically discuss the implications of our 
analysis. 

PACS numbers: 04.50.Kd, 95.36.-l-x, 98.80.-k, 98.80.Cq,11.25.-w 


I. INTRODUCTION 


Bouncing cosmology provides us with an appealing solution of the initial singularity problem, which is a 

rather unwanted feature in cosmological theories. This is because, in the context of bouncing cosmology, the Universe 
contracts until a minimal radius is reached, and after that point it expands. Therefore, the Universe never collapses to 
a singular point, thus avoiding the initial singularity. The cosmological bounces can appear in two main categories, the 
Loop Quantum Cosmology 111 matter bounce theories [l3 and also in theories that make use of scalar fields in order 
to avoid singularities and generate bounces [MIS- In addition to these, modified gravity also offers a consistent 
description of bouncing cosmology fl^.[l3l|. In the context of bouncing cosmology, the acceleration and thermal history 
of our Universe are consistently described (see Refs, for review on this issue) and also certain CMB anomalies 

on large angular scale find a satisfactory explanation |15l | . In order for a bounce to occur in standard Einstein-Hilbert 
gravity, the null energy condition has to be violated, something that is only possible for physical systems for which the 
Hamiltonian is bounded from below [T^ . Some drawbacks that come along with the appealing features of the bounces 
are that, bouncing cosmologies suffer from ghost and primordial instabilities during the contracting phase, rendering 
the contracting phase a problematic era during the evolution. The issue of the contracting phase is known as BKL 
instability 0. Ghost instabilities can be resolved in the context of Galileon and ghost condensate models Sd, 
and the issue of the contracting phase is successfully resolved in the context of the ekpyrotic contraction theories 
Bid- Actually, as was shown in Ref. [l^ , the BKL instability issue can be successfully resolved in the context of 
an ekpyrotic contraction. 

On the other hand, the initial singularity is not the only type of singularity that may occur in a cosmological 
theory, since there exist other types of milder singularities, which are called sudden [20l - [^ or finite time singularities 
[ 2 ^ . The initial singularity is a crushing type singularity, in which the strong energy theorems apply and 

geo desics incompleteness occurs at these points. However, in the context of finite time singularities, only the Big Rip 
is of crushing type, with the rest of the finite time singularities being milder, and also geodesics incompleteness 
does not necessarily occur for the non-crushing types singularities. In the context of general relativity, crushing type 
singularities occur during the process of gravitational collapse, and several conjectures have been proposed that point 
out the need to protect the rest of the Universe from these naked singularities with a “cloth”. This is the cosmic 
censorship hypothesis [30l| - which up to date has not been proved yet. For an informative account on this issue see 
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[3ll - l^ . In cosmology however there is no way to “dress” these singular points, therefore the complete understanding 
of their nature and implications is compelling. It seems that the non-crushing type singularities offer a good testing 
ground for the complete understanding of finite time singularities and their consequences in the cosmological evolution 
of our Universe. In Refs. [33 - 1^ . we studied the Type IV finite time singularities and their implications in the 
cosmological evolution (especially, after inflation), in the context of single scalar [s^, or multiple scalar fields 
In this pap er we shall study how a Type IV singular bouncing cosmology can be generated from a general pure F{R) 
[14 . I37h 41| gravity, with pure meaning that no matter fluids are assumed to be present. An important assumption 
we shall make is that the bouncing point coincides with the point that the Type IV singularity occurs. We shall be 
particularly interested in finding the F{R) theory which generates the bonnce near the bouncing point, which is also 
the point at which the Type IV sing ularity occurs. For a list of review and important papers on F{R) theories of 
gravity, the reader is referred to [IJ, j^. In order to reveal which Jordan frame F{R) gravity can successfully 
describe the singular bounce near the Type IV singularity, we shall use some very well known reconstruction techniques 
[33 - l4lj| . and as we demonstrate, the resulting Jordan frame F{R) gravity is an gravity plus cosmological constant. 
Having found this form of the F{R) gravity, we also investigate the Einstein frame implications of this F{R) gravity, 
assuming a quasi de Sitter solution in the Einstein frame. The observational implications of the Einstein frame 
canonical scalar theory are also studied in detail. Finally, we also study the behavior of the effective equation of state 
(EoS) corresponding to the Type IV singular bounce solution and compare the resulting picture with other bouncing 
cosmologies. 

This paper is organized as follows: In section II we present in brief all the essential information with regards to 
finite time singularities and also we discuss our motivation to use a Type IV singularity for our study. In section 
III we study the general properties of the Type IV singular bounce we shall use in the forthcoming sections and in 
section IV, using well known reconstruction techniques, we investigate which F(R) gravity can successfully describe 
the Type IV singularity near the singularity point, which we chose to coincide with the bouncing point. Also the 
stability of the resulting solution is discussed in the end of the section. In section V, we study the Einstein frame 
canonical scalar-tensor theory corresponding to the Jordan frame F{R) gravity we found in section IV. In addition, 
we also investigate the observational indices of the Einstein frame scalar theory, assuming a quasi de Sitter solution 
or a de Sitter solution exists in the Einstein frame. In section VI, we perform a full analysis of the EoS corresponding 
to the Type IV singular bounce under study, and also we compare the behavior of the Type IV singular bounce, to 
other bouncing cosmologies. Finally, a critical discussion on the results along with the concluding remarks follow in 
section VII. 


Geometric background conventions 

Before we start, it is worth mentioning the background geometric conventions we shall use in this article. We 
initially work in the Jordan frame for all F(R) gravities we shall discuss, and we adopt the metric formalism approach 
[j^l . In addition, we assume that the background geometry consists of a pseudo-Riemannian manifold, which locally 
is a Lorentz metric, a flat FRW one, with line element, 

ds^ =—dt"^ + ( 1 ) 

i 

In this geometric background, the Ricci scalar reads, 

R = 6(2ij2 + H), (2) 

with H{t) denoting as usual the Hubble rate, and the “dot” denotes differentiation with respect to the cosmic time 
t. Finally, we choose the affine connection on this manifold to be the Levi-Civita one, which is a metric compatible, 
symmetric and torsion-less. 


II. FINITE-TIME SINGULARITIES ESSENTIALS 

The classification of finite-time cosmological singularities was extensively done in a formal way in Refs. [13, 
and we briefly recall the basic features of this classification, adopting the notation of Refs. [13, llg|. There are four 
types of finite time cosmological singularities, the Type I,H,HI and Type IV singularities, and these are classified in 
the following way[2l[2g, 

• Type I (“Known as Big Rip Singularity”) : This type of cosmological singularity is the most severe among 
finite time cosmological singularities and it is a singularity of crushing type. It occurs when as the cosmic time 
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approaches a specific time (tg), that is when t —?> tg, the effective energy density peff, the scale factor o(t), and 
also the effective pressure Peff diverge, that is, a —>■ oo, peff oo, and |peff| oo. For an important stream of 
papers with regards to the Big Rip singularity, the reader is referred to. Ref. [ 2 ^ [H, [1^ 

• Type II (Known as “Sudden Singularity”) [2l|, [ 2 ^ : This singularity occurs when, as the cosmic time approaches 
t —>■ tg, only the scale factor a and the effective energy density peff take bounded values, that is, a —>■ Ug, 
PeS —>■ Ps, with both as,Ps < 00 , with the effective pressure diverging as t — >■ tg, that is, |peff| 00 . This case 
for example occurs when the second and higher derivatives of the scale factor diverge. 

• Type III : This singularity occurs when, as the cosmic time approaches t —>■ tg, only the scale factor remains 
finite a —>■ Og, but both the effective energy density and the corresponding effective pressure diverge, that is, 
|peff| 00 and peff —>■ 00 . Equivalently, this means that the first and higher derivatives of the scale factor 
diverge. 

• Type IV : This type of singularity is the most mild among all the three aforementioned types of finite time 
singularities, and we shall focus on this type of singularity in the following. For a detailed study on this 
finite time singularity, see Ref. [ 2 ^. This singularity occurs when, as the cosmic time approaches tg, all the 
cosmological physical quantities remain finite, that is, the effective energy density Peff —>■ Ps, the effective pressure 
|peff| —>■ Ps and the scale factor a —>■ Og, but the higher derivatives of the Hubble rate diverge. 

In the next section we shall further analyze in brief our motivation to use the Type IV kind of finite time cosmological 
singularity for our analysis of the singular bounce. 


A. Why choosing a Type IV singularity: Brief discnssion 

As we already mentioned, among all the types of finite singularities which we presented previously, the most mild is 
the Type IV, with mild referring to the geodesics incompleteness issue that might occur at finite time singularities. The 
crushing type singularities, such as the Big Rip or the initial singularity (T7j |. always lead to severe phenomena, 
which when considered classically, these are flaws of the theory, like the singularity in the Coulomb potential in 
classical electrodynamics. For these singularities, the energy theorems are violated, therefore these can be seen as 
either flaws of the theory or indicators that a new, quantum maybe, theory describes the physical system at the 
energies that these correspond. 

On the other hand, the Type IV singularity is less harmful since no geodesic incompleteness necessarily occurs and 
also all the energy theorems of Hawking and Penrose are satisfied, so these singularities are not so severe. However, 
their presence and possible consequences should be scrutinized in order to fully understand what these singularities 
indicate. In this paper we shall study a bouncing cosmology that has in its evolution a Type IV singularity. Notice 
that in bouncing cosmology no initial singularity appears, so the only singularity that occurs during the bounce 
evolution is the Type IV. These appealing properties of the Type IV singularity motivated us to realize a non-singular 
bouncing cosmology by using a Jordan frame F{R) gravity. 


III. DETAILED DESCRIPTION OE THE BOUNCE SOLUTION 

In this section we shall present the Type IV singular bouncing cosmology that we shall extensively study in the 
following sections. Before we start, it is worth recalling the basic properties of the bouncing cosmologies. For details 
on this see Refs. M- 

A cosmological bounce consists from two eras, a contraction era and an expansion era. At the beginning, and 
during the contraction, the scale factor decreases, that is d < 0, until the Universe reaches a minimal radius, where 
d = 0, it bounces off and starts to expand, with d > 0. The fact that the Universe reaches a minimal radius is what 
renders the bouncing cosmology so appealing, since the initial sin^larity is avoided in this way and there exists also 
the possibility of describing successfully early time acceleration (JH, [13 > thus avoiding the standard description with 
inflationary models. 

Let us assume that the bouncing point occurs at a cosmic time U. When the Hubble rate is taken into account, in 
the contracting phase, which means that t <ta the Hubble rate is negative H{t) < 0, at the bouncing point becomes 
equal to zero , H{ts) =0, while for t > tg, the Hubble rate is positive, H{t) > 0. 

The singular bounce we shall consider in this paper, has the following scale factor, 


a{t) = e- 




( 3 ) 
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where we normalized the scale factor to be equal to one at the bouncing point, that is a(ts) = 1, and also e and /o are 
constant parameters, with the values that e is allowed to take to be specified shortly. The Hubble rate corresponding 
to the scale factor ([3|) is equal to. 


= (4) 

Note that in the Planck unit system, the Hubble rate is measured in eV, the time is measured in (eV)“^, so the 
parameter /o is measured in However, we shall express time in seconds, so the Hubble rate is measured in 

(sec)“^ and consequently the parameter fo is measured in (sec)“^®“^. We adopt these units in the rest of this paper. 

Taking into account the classification of singularities we presented in the previous section, the Type IV singularity 
occurs when the exponent in Eq. o, is 2£: +1 > 1, which means that e > 0, so for all positive values of the parameter 
e. However, we shall assume that e < 1 in which case, the bounce ([3]) is a small deformation of the well known [l^ 
bouncing cosmology. 


a(t) 


(5) 


Notice that the Type IV singularity occurs at t = tg, which is exactly the bouncing point. 

It is worth elaborating on the Type IV singularity caused by the structure of the Hubble rate (jl]) . Set for simplicity 
2e + 1 = /3, so that the Hubble rate becomes, 

H{t) = 2(1 + e)fo {t - tsf . (6) 


Let us see for which values of /3 the Type IV occurs. In the list below we quote all the possibilities for a finite 
singularity to occur, for various values of the parameter /3, 

• /3 < —1 corresponds to the Type I singularity. 

• — I < /3 < 0 corresponds to Type HI singularity. 

• 0 < /3 < 1 corresponds to Type H singularity. 

• /3 > 1 corresponds to Type IV singularity. 

The Type IV singularity case occurs when /3 > 1. This can be easily seen, since in the Type IV case, the singularity 
occurs when the higher derivatives of the Hubble rate are divergent, which means that, 


d”iJ(t) 

dt” 


—>■ oo, 


(7) 


for some n > 2. Let us compute the lowest derivative for which the Type IV singularity could occur, which is for 

n = 2, 


^^=2(I + £)/o/3(/I-I)(^-^g)''-^ (8) 

which is clearly divergent when l</3<2. If/3>2, then the second derivative of the Hubble rate is finite, but the 
third derivative becomes divergent, 

= 2(1 + e)fo/3{P - l)(/3 - 2) (f - tsf-^ . (9) 

In order for the bounce (HI) to be a deformation of the bounce ([5|), we assume that e < 1. Since ,5 = 2e + 1, this means 
that certainly /3 > 1 but also that /3<2, sol</3<2, which means that the second derivative of the Hubble rate 
is divergent (see Eq. ([5])). Consequently a Type IV singularity occurs for 0 < e < ^ and for these values of e, the 
bounce (l3|) is a deformation of the bounce ([5]), since e < 1. We have to note that we also restrict ourselves to positive 
values of the parameter e. 

Before we proceed, we need to address another issue where a possible inconsistency might occur. Throughout the 
article we shall assume that 0 < £ < | and the values that e is allowed to take are chosen is such a way so that the 
scale factor and the Hubble rate never become complex. With regards to the scale factor, this would require that the 
exponent oi t — tg in Eq. takes the following form, 

2{s+l) ={t- , 


{t - tg) 


( 10 ) 
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with n and m being arbitrary integers appropriately chosen so that 0 < e < ^. One convenient choice that never 
makes the scale factor complex but also renders the parameter e smaller than one, is for n = 12 and m = 5. In this 
case, the scale factor (jS]) reads. 


which is never complex, for all For the choice of n and m, the parameter e reads e 
hereafter. For e = ^, the Hubble rate reads. 


( 11 ) 


jj, so we use this value 


H{t) = 2(1 + e)fo {t - = 2(1 + e)/o {t - {t - ts) = {t - ts) = 2(1 + e)fo {t - ts) {t - 4) , (12) 

which clearly never takes complex values (for the same reason as in the scale factor case, since {t — tg)^ = 
1 

((t — “), but can become negative for t < tg due to the last term {t — tg). Having cleared out this vague 

spot, let us see now how the scale factor and the Hubble rate behave as functions of cosmic time. To this end, we 
choose for illustrative purposes only, tg = 10“^®sec, e = and /o = 0.001(sec)“^®“^. In Fig. [H we plot the time 




FIG. 1: The scale factor a{t) (left plot) and the Hubble rate (right plot) as a function of the cosmic time t, for ts = 10 ®®sec, 
£ = i and /o = 0.001(sec)-2=-^ for a{t) = 

dependence of both the scale factor ([3]) and of the Hubble rate (j4]). As we can see all the qualitative features of the 
bounce are satisfied, that is, before the bounce H{t) < 0, at the bounce H{t) = 0, and after the bounce H{t) > 0. 
In addition, the contraction and expansion can be observed by looking the behavior of the scale factor, since for 
t < tg, the scale factor decreases until a minimal value is reached at t = tg, and then for t > tg the Universe expands. 
Finally, in Fig. O we plot the behavior of the second derivative of the Hubble rate H"{t) as a function of time, for 



FIG. 2: The second derivative of the Hubble rate H"{t) as a function of the cosmic time t, for ts = 10 ^®sec, £ = ^ and 
fo = 0.001(sec)“^'^“^, for the bounce a{t) = ^ ^ ' 

tg = 10“^®sec, £ = and fo = 0.001(sec)“^'^“^. As we can see, at the Type IV singularity, the function H"{t) 
blows-up for t ^ tg, as expected. 


1 


Note that if we take the scale factor to be equal to a{t) 


e \ 


24 

, then for t < tg, the expression of the scale factor contains 


(—1)TT^ which has one negative but real branch and two complex branches, that is (—1)TT = —1,0.959493 + 0.2817332,0.959493 — 
0.2817332, so by keeping the real branch, we obtain again a real number, since the resulting scale factor contains (—1)^"^, which is 
positive 
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IV. BOUNCE SOLUTION FROM JORDAN FRAME F{R) GRAVITY 

The F{R) modified gravity theoretical framework can realize cosmological scenarios, that were, to some extent, 
exotic for ordinary general relativity. Noteworthy is the fact that reconstructed F(R) gravity models can have a 
direct relation with viable gravity models based on the Khoury chameleon scenario [42j| . In this section, by using well 
known reconstruction techniques [sol - l^ . we shall investigate which F{R) gravity can generate the Type IV singular 
bouncing cosmology of Eq. The focus will be given for times near the Type IV singularity, that is when t —>■ 
For important reviews and papers on the reconstruction techniques we shall use, consult [37|,l38| and references therein. 
Consider the following Jordan frame F{R) gravity, with action, 

5 = ^ / d^Xy/^F{R) + Sm , (13) 

where k, is related to Newton’s constant = SttG, and in addition Sm denotes the action of all matter fluids present. 
In the context of the metric formalism of F{R) gravity (3^ . upon variation with respect to the metric we obtain 
the following equations of motion, 

^ + V^V,F'{R) - g^,UF'{R)yj . (14) 

where the prime denotes differentiation with respect to the curvature scalar R and as usual, T^j/ is the energy 
momentum tensor that receives contributions from all matter fluids. By observing Eq. da, it is obvious that the 
energy momentum tensor receives an extra contribution which originates from the F{R) gravitational sector. This 
extra contribution, sometimes called geometric contribution, is what makes F{R) gravity a modified theory of standard 
Einstein-Hilbert gravity. Particularly, the extra contribution to the energy momentum tensor is equal to, 

T^eff _ 1 f F{R) — RF'{R) ^ ^ ^ ^ 

Rfiu ~ ^2 ( 2 ^fi^vF (i?) gfj^i/DF (i?) 1 . (15) 

In addition, we assume that the spacetime metric is the flat FRW of Eq. ©• 

What we are interested in now, is to find which pure F{R) gravity can generate the cosmological evolution of the 
singular bounce with scale factor given in Eq. ©■ With pure F{R) it is meant that no matter fluids are assumed 
to be present. We make use of quite well known reconstruction techniques [39l - l4ll| and we shall focus our study for 
cosmological times near the Type IV singularity, which occurs ai t = tg. Start from a general pure F{R) gravity in 
the Jordan frame with action, 

5 = ^ / d^x^F{R). (16) 

The first FRW equation can be obtained upon variation of action (1161) . with respect to the metric tensor and 
reads, 

- 18 (4H{tfH{t) + F"{R) + 3 + H{t)) F'{R) - = 0 . (17) 

The reconstruction technique we shall use, involves an auxiliary scalar field </>, which enters the action of Eq. m, in 
the following way, 

S = J d\./^ {P{^)R + Qm . (18) 

In the absence of a kinetic term for the scalar field, this is practically a non-dynamical degree of freedom, this is 
why we called it an auxiliary degree of freedom. By varying the action of Eq. m, with respect to (j), we obtain the 
following algebraic relation. 


P'(<^)R + Q'(<).)=0, (19) 

with P'(<^) = dP{<j))/d(j) and Q'{4>) = dQ{(j))/d(j). This equation is of particular importance, since it will provide us 
with the function 4>{R)., if it can be solved explicitly with respect to the scalar field. This is a very important point 
for the reconstruction method, since by finding (^(R), we can substitute it in the following relation, 


F{cl){R))=P{^{Rj)R + Q{^{R)). 


( 20 ) 
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and have explicitly the F{R) gravity. So what we need to know is the explicit form of P{(f>) and Q{(f>) and in order 
to find explicitly their functional form, the action (1181) must be varied with respect to the metric, so we obtain the 
following relation, 

0 = - GH^Pit) - Q{t) - = 0 , 

0 = (4H + 677^) Pit) + Qit) + = 0 . (21) 

Then, by eliminating the function Qi4){t)) from Eq. (HU, we get the following second order differential equation, 

2^^®-2i7(t)^:^+4i7P(t) =0. (22) 

at^ at 

Notice that in Eqs. (HU and (l2^ . we used the cosmic time t as a variable of the functions P(^!>) and Q{(l)), instead of 
(j). Practically, these variables are identified within the reconstruction technique we are using, and the relation cj) = t 
is valid for a wide range of field values, since the actions of Eqs. (1161) and (I18|) . are equivalent from a mathematical 
point of view. Eor details on this account, see the Appendix of Ref. (39j. So, for a known cosmological evolution with 
a specified Hubble rate, the solution of the differential equation Eq. (HU, yields the function P{t), and by substituting 
the result to Eq. HU, we obtain the function Q{t). 

Let us proceed to find an analytic form for the function P{t) and then proceed to find the E(i?) gravity near the 
Type IV singularity. By using the Hubble rate (lU, the differential equation (1^ reads, 

2— ^ — 4/o(l + e) (t — _|_ 8/o(l + £)(1 + 2e) {t — tg)^^ P{t) = 0. (23) 

Since we are interested in finding the E(i?) gravity near the Type IV singularity, it is worth changing the variable t, 
to the new one x, defined to be a: = t — tg. In this way, as t approaches the singularity, the variable x approaches 
zero, or schematically as t —>■ tg, then a; —?> 0. Using x as a variable, the differential equation (1251) becomes. 


^ ^ 2 ^^ - 2/o(l + e)x^^+^+ e)(l + 2e)x^^P(x) = 0. (24) 

which can be analytically solved by setting z = and hence it becomes, 

(2e + 2) z — ^^2 + (2s + 2)(— 2/o(l + e)z + 2e + 1) —+ 4/o(l + e)(l + 2e)P(z) = 0, (25) 

The solution to this differential equation is, 

P(z) = (2z + 2ze)TOCi[/(^-^i^,l + ^^^,/oz) +(2z + 2ze)^C2L;r(/o^) , (26) 


where the functions U{a,b,z) and F^iz) are the confluent Hypergeometric function and the generalized Laguerre 
polynomial respectively, the parameters Ci, C 2 are constant arbitrary real numbers and finally m and n stand for. 


l + 4e 1 

— _ 777 ^ _ 

2(1 + £)’ 2(1+ £) 


(27) 


In terms of the variable x, the function P{x) reads. 


P(x) = (2x2'^+2 + 2x2"+2£)mrc'^[/ 
Then, the function Q{x) easily follows. 


1 + 4£ 

2(1+ s) 


> 1 + 


1 


Qix) =3 2^+2 + 2e/oX^+^'^(l + £) (l+£) 2 + 2 e X 


2(1 + £) 

■ /o(1 + 4£)C: 


,/ox2^+^ ) +(2 x2''+2 + 2x2^+^£)51t+^C'2L™ ■ 

(28) 


1 + £ 


-K 


1 + £ 


- 2 ( 1 +U ( - (1 + 4/ox3+^^(1 + £)2))Ci) p( 


l-2s , 
2 + 2 £’‘ 
1 + 4£ 


1 


2 + 2 £ 




2(l+e) + l 


2 + 2 £ 


, 1 + 


1 


2 + 2 £ 




,2(l+e) + l \ 


+ C 2 (2/ox2+2-(l + £)L-^ - (1 + 4/ox3+4-(l + e)^) 


( 29 ) 































where we introduced the variables ni^mi which are equal to, 


and also 712 and m 2 stand for, 


-1 + 2 £ ^ 1 

= 5(TT7)' ”« = i + 2T 


712 


1 +4£ 

2 + 2e ’ 


m 2 


1 

2 + 2e ■ 


( 30 ) 


(31) 


Having at hand the functions P(x) and Q(x), in principle by substituting in Eq. (IT^ . we can obtain the functional 
dependence of x as a function of R. However, the form of the functions P{x) and Q{x) does not allow an easy analytic 
manipulation of the resulting algebraic equation. However, since we are interested in the behavior of the F{R) near 
the bounce, and since this limiting case is reached in the limit x —>■ 0, we shall approximate P'{x) and Q'(x) near the 
bouncing point. By taking the derivative of Eq. (I28L approximating the resulting expression in the limit x —>■ 0 and 
by keeping only the leading order terms, the function P'{x) reads, 


P'(x) 


A 

^2(e+l) + l 


0(x), 


(32) 


where we introduced the constant parameter A, which can be found in the Appendix B. 
By doing the same to the other function Q'{x), we obtain in the small x limit, 


«'W-j5(;li)Tr+c + oW. (33) 

where the parameters B and C can be also found in the Appendix B, and also notice that only leading order terms 
where kept. Upon substitution of Eqs. (1321) and (l33t into Eq. (HU), we obtain the function x{R), 


X ~ - 


C 


AR + B 


2{e+l) + l 


(34) 


Having x{R) at hand, the F{R) gravity near the Type IV singularity easily follows by making use of Eq. (1^01) . so the 
resulting expression of the F{R) gravity is. 


BA 

F(R)c,-^R^-2^R--+C. 


(35) 


Therefore, the F(R) gravity that generates the bounce (|31) near the Type IV singularity is a nearly R^ gravity . We 
can bring the resulting expression to be exactly an Einstein-Hilbert gravity plus curvature corrections by appropriately 
choosing the free parameter Ci to satisfy the following constraint, 


which holds true if Ci is chosen to be. 


Cl =- 


2^ = 1. 
C 


(1 + 2(£ + i))2r(^t|g^)r( wg)) 

I2/o(I + 32(£+I)) 


(36) 


(37) 


((2 + 2(e + i))r(^tw)r(i + ,^5^) - 2(1 + 22(£ + i))r(^tw)r(2 + 

The resulting F{R) gravity near the Type IV singularity then reads. 


£+ 1 ) 


42 

F{R)c^R-—R^- — +C. 


■ l + 2(e + l) 


(38) 


Recall that we assumed e <C I, and by also assuming that /o > 0, it can easily be shown that the parameter Ci < 0. 
Consequently, the parameter C, defined in Eq. (II04L is negative, therefore the coefficient in front of the R^ in Eq. 
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(l38l) is positive. It is worth redefining the coefficients of the F{R) gravity, for later convenience, so we introduce the 
following new variables, 

C 

^0 = -^^, A = +C . (39) 

Notice that, since C < 0, the coefficient Cq is positive. In terms of the new coefficients, the Jordan frame F{R) gravity 
for the bounce ([3]), near the Type IV singularity has the following form, 

F{R)=R+^^+A. (40) 

This Jordan frame F(R) gravity has a particularly interesting Einstein frame counterpart, since it corresponds to an 
nearly R^ scalar theory in the Einstein frame. We shall discuss this issue in detail in a later section. 


A. Stability analysis of F{R) gravity solution 


Having at hand the reconstructed F{R) gravity near the bounce, it is a straightforward task to examine the 
stability of the solution near the bounce. What is expected is that the system of differential equations, when viewed 
as a dynamical system, is unstable near the bounce, since the bounce is not an ending state of the system but just a 
passing point during the Universe’s evolution. 

In order to address formally the instability issue, we shall use the approach adopted in Ref. [39l - l4lj |. The starting 
point of our analysis is the following equation. 




dF ^ dt 

where in our case, = {(j) — The expression in Eq. (HD) can be written as follows, 

1^ 

dF . 


(41) 




dcj) dP 

and accordingly can be recast as follows, 


^ + m] ^ 0 . ( 43 , 


d02 (^(j) 

We introduce the function S to be equal to, 


4 = 

dt 


(44) 


Practically, the parameter 6 measures the exact way that perturbations behave for the solutions we presented in the 
previous section, with regards to P(0), if the system of equations are treated as a dynamical system. In terms of S, 
Eq. (H51) can be written as follows. 




with u}(t) standing for. 


;(t) = 2 




™ + 2g'(,)))P(0) 


(f)—t 


(45) 


(46) 


Thereby, if (^ > 0, the dynamical system of Eq. (HSl) is stable, since the coefficient of the dynamical variable 5 is 
rendered negative. However, if w < 0 the dynamical system is unstable, since the perturbations grow in an exponential¬ 
like way. Since we are interested in the limit t ^ tg, the functional form of uj{4i), near the Type IV singularity, where 
the bounce occurs, reads. 


;(</,) ~ 


2 (6 + £ (7 - ((0 - tgfY {2 + x-de) + 2e)) 
((/.-U)(-3+(-2 + 4((</>-U)2)^)£) 


(47) 
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which for 


ts can be further approximated by the following expression, 




2(6 + £(7 + 2£)) 
{(j) - ts){-3 - 2e)' 


(48) 


As is obvious by looking at Eq. (l4^ . w((/)) is negative for cj) > ts, while it is positive for (j) < ts, so the system is 
conditionally unstable (as we anticipated), since the point (f> = tg is a saddle point of the dynamical system (1451) . 


V. EINSTEIN FRAME ANALYSIS OF THE F{R) BOUNCE 

In the previous section we demonstrated that the Jordan frame F{R) gravity of Eq. (HOI) can realize the bounce 
cosmology of Eq. ([3]). In this section we show that the Einstein frame counterpart of the F(R) gravity (l40l) . is a 
deformed form of R^ inflation in the Einstein frame (for a similar approach to ours, see Ref. [45l|l. Particularly, it 
corresponds to a canonical scalar theory with a nearly R^ potential (43| . This kind of potentials were first studied in 
( 4 ^ . We have to note that we do not assume that we start in the Jordan frame with the scale factor of Eq. (|31), but 
we start with a convenient metric that, when conformally transformed to the Einstein frame, it can generate nearly 
Starobinsky inflation. We start off with the action of the F{R) gravity (HU)) . 

+ ^ + ( 49 ) 

where 5 ^ 1 / is the Jordan frame metric tensor. We introduce the auxiliary field A, so in terms of this auxiliary scalar, 
the action can be written as follows, 

5 = ^ / d^xv^ {F'{A){R -A)+ F{A)) . (50) 

Upon variation with respect to the scalar A, the solution A = i? is obtained, thus the mathematical equivalence of 
the actions and (|50l) is verified. By using the canonical transformation. 


If = 


2k2 


ln(F'(A)), 


(51) 


we are transferred to the Einstein frame. In Eq. dSB, the scalar field (f is the Einstein frame canonical scalar field. 
Conformally transforming the Jordan frame metric 5 ^ 1 ,, 


g^iv — e ^dfiv ■, 


we obtain the Einstein frame scalar field action, 
S = 




A F(A) 
F'(A) ~ F'(Ay 


with the scalar potential V{(p) being equal to 


-y/2K^/3(p'^ _ p p ^g-\/2K^/393^ ^ 


F'{A) F'{AY 2k2 

For the Jordan frame F{R) gravity (^(11) . the canonical scalar field potential of Eq. (I54L becomes 

V{ip) ~ Co + + Cie-v^'"'^ , 

where the constant parameters Ci and C 2 are related to Cq and A that appear in Eq. (HOI) as follows 


(52) 


(53) 


(54) 


(55) 


Ci=-2Co, C 2 = Co-A. 


( 56 ) 
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The canonical scalar field potential (l55l) corresponds to a nearly E? inflationary potential, which in the case that 
the coefficient C 2 is equal to Cq, then the potential becomes exactly the model potential in the Einstein frame, 
that is. 


= (57) 

This by any means does not imply that the bounce in the Jordan frame with metric that has a scale factor ©, 
corresponds to an inflationary solution when viewed in the Einstein frame. That would require a very special conformal 
transformation that may not necessarily lead to a de Sitter or even quasi de Sitter solution in the Einstein frame. An 
interesting scenario occurs when the Jordan frame is such that, when conformally transformed according to the 
transformation (EH), it becomes a de Sitter or quasi de Sitter metric in the Einstein frame, with the scalar potential 
being that of Eq. (j55|) . This is a quite appealing scenario and in the next section we shall consider the observational 
indices of this nearly model (1551) . Some interesting studies related to conformal transformation between frames 
and singularities can be found in [44| . 


A. Observational indices of the Einstein frame F{R) gravity 

In this section we shall study the observational implications of the Einstein frame canonical scalar field model with 
potential that of Eq. dSSl). The potential of Eq. (l55l) in view of the constraints (IMl) . reads, 

V{^) = Co + Cse-^v^'^^ - 2Coe-'/l’'=‘^. (58) 

For the canonical scalar theory in the Einstein frame, assuming a flat FRW metric, the energy density and pressure 
of the model are equal to, 


• 2 *2 

Pv = \ + = Y - ^(t’), (59) 

where as usual, the ’’dot” indicates differentiation with respect to the cosmic time. The Friedmann equations in the 
presence of the canonical scalar </?, are given by, 

^ = Y + n7’), (60) 

and also the following second order equation of motion for the scalar field ip is satisfied, 

(/3 + =—E'((/?). (61) 


Before getting to the detailed calculation of the observational indices for the model (1551) . it is worth presenting in 
brief the essentials of the slow-roll app roximation in inflationary theories. For a detailed presentation of these issues, 
the reader is referred to Refs. [11,143. Historically, the slow-roll approximation [4l,|4l|, was introduced in order to 
solve the graceful exit of inflation problem and nowadays is frequently used in most inflation predicting theories. The 
slow-roll conditions are based on the following constraint for the canonical scalar field, 

« V {^), (62) 

where it is assumed that this constraint holds true for a sufficiently long period of time, with the latter feature being 
model dependent nevertheless. The constraint of Eq. (|62|l is called the first slow-roll condition, and it guarantees a 
long and finite acceleration era. Furthermore, another constraint has to be imposed on the first slow-roll condition, 
so that the slow-roll accelerating era lasts for a sufficiently long period of time, which is the following, 


|v3| < 


dV{v) 

dip 


(63) 


The condition of Eq. ((63ll is called the second slow-roll condition, which by taking into account the equation of motion 
of the canonical scalar field, namely Eq. (ED), in a flat FRW background, it becomes as follows, 


\(p\ < 37J|(yj|. 


( 64 ) 
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In view of the two slow-roll conditions, namely Eqs. (EH) and (l64l) . the equation of motion of the scalar field (IhTI) 
becomes, 


1 dV{ip) 

dip 


(65) 


while the FRW equations of Eq. (pl) read, 

~ ZHif ~ -V'{ip). 


( 66 ) 


The two slow-roll conditions, namely Eqs. (l63ll and (l64ll . can be shown that they are equivalent to the following two 
relations, 


v'iv) 

V{^) 


< 2k^, 


which accordingly can be recast as follows. 



< , 


e < 1, < 1, 


(67) 


( 68 ) 


These two conditions (IS51) constitute the slow-roll conditions. The parameters e and rj are known as the slow-roll 
parameters, and these, in the context of slow-roll approximation, are formally defined to be equal to, 


V{(^) ) ’ / 


(69) 


Notice that the prime in all the equations above denotes differentiation with respect to the canonical scalar field ip. 
The observational indices that are currently scrutinized by the Planck collaboration (50j |. are written in terms of the 
slow-roll parameters. Particularly we shall be interested in the spectral index of the primordial curvature fluctuations, 
denoted as rig, and the tensor-to-scalar ratio, denoted as r, which are expressed in terms of the slow-roll parameters 
(IMD in the following way , 


rig ^ 1 — 6 e -I- 277 , r = 16e. (70) 

Finally, the latest Planck (2015) observational data predict for the observational indices of Eq. (I70|) . the following 
values. 


rig = 0.9655 ±0.0062, r < 0.11. 


(71) 


After this brief introduction to the slow-roll approximation and related observational indices, we now proceed to the 
phenomenological analysis of the Einstein frame canonical scalar field theory counterpart of the F(R) gravity (I40L 
with Einstein frame canonical scalar potential (1551) . The model of Eq. (1551) was also presented in Ref. [^, to which 
reference we refer the reader for further details. As it can be shown, the minimum of the canonical scalar potential 
(1551) . is at = 0, since the critical point of the equation V'{p) = 0 is actually <p = 0. Also, since the second derivative 
of the potential at (/? = 0 is. 


E"(0) = 


3 ^ 3 


(72) 


which is positive when C 2 > the critical point (/? = 0 is a global minimum of the potential V{p). We recall here 
that the inflationary evolution of the canonical scalar field p goes as follows: during the inflationary era, the scalar 
field has very large values {p) and inflation ends when the slow-roll parameters of Eq. (I69|) . become of the order one. 
Eventually, the field tends to the value p, which corresponds to the minimum of the potential. For the potential (j58L 
the Eqs. (1651) and (l 66 l) in the slow-roll approximation become. 


3772 ^ _ 2e-V^'=‘^7 ± 4C'2e-2'/f«<^K2 


, mp^- 


- 4C'2K^) 


4k2 ’ ^ 

where for convenience we introduced the parameter 7 , which is defined in terms of Co as follows, 

7 = 4«;^Co. 


(73) 


(74) 
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Since during inflation, the canonical scalar field has large values {ip —>■ oo), the FRW equations and the scalar field 
equations (El become, 


where we ignored subdominant terms and more importantly we assumed that, 

7^ » (76) 

which ensures that a quasi de Sitter solution can be achieved. So during inflation, where ip ^ oo, the slow-roll indices 
read, 


e 



rj ~ 


4 




(77) 


The inflationary era ends for the value of p, for which the parameter e (or rj) becomes of order one, that is, e{pend) — 1, 
which for the parameter e appearing in Eq. El, occurs for, 


Pend — 


3 1 

2 K 


In 2- 


(78) 


By using Eq. (El, the observational indices can be written as follows. 


n, ~ 1 




64 


6 - 3evT'=v ^_2 j ^ 3 ^ 


(79) 


3 -2-heV3' 


We can write these by using the e—folding number as follows (see also 0):, 

2 12 

In order to see this, recall that the e—folding number in the slow-roll approximation is defined to be equal to. 


N 


V{^) 

V'{p) 


Ap, 


(80) 


(81) 


where p*, is an initial value of the scalar field p which is assumed to be p* ^ Pend- By using this approximation, we 
obtain that. 


N 


3 /T, 

= -eV 3‘ 


if K 


(82) 


for which value, the slow-roll parameters become. 


e 



rj ~ 




(83) 


and by combining Eqs. (IMI) and (IMl) . we easily obtain Eq. (1^ . Then, by looking Eq. (|M)l . in order to achieve for 
exarnple iV = 60 e—folding, it is required that the initial value of the scalar field is approximately, pi ~ (see 
also [ 4 ^ for details). The resulting picture is compatible with the 2015 Planck data [H^, since for N = 60, we obtain 
that ns ~ 0.9665, and r ~ 0.0029. In conclusion, as was also pointed out in detail in 0|, the resulting picture is 
pretty much alike the standard model, with the only difference being traced on the fact that the minimum of the 
potential is non-zero. Therefore, the F{R) gravity model of (HUl) . which describes a Jordan frame R^ gravity plus 
cosmological constant, in the Jordan frame can give rise to the singular bounce of Eq. (0, and when the metric in 
the Jordan frame is appropriately chosen, the same model can give rise to an Einstein frame inflationary potential 
which is a modification of the R^ model, with observational indices compatible with current observational data. 
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VI. THE EFFECTIVE EQUATION OF STATE FOR THE SINGULAR BOUNCE AND COMPARISON 

WITH NON SINGULAR BOUNCE 


Having studied the F{R) gravity that can generate the Type IV singularity, in this section we shall study in detail 
the EoS corresponding to the singular bounce of Eq. ©■ Eor both a standard Einstein-Hilbert gravity background 
and a modified gravity background^, the EoS for a Hubble rate H{t) is defined to be [3^ . 


= -1 - 


2H{t) 


so for the Hubble rate of Eq. 0, the EoS reads, 


WeS = -1 - 


(1 + 2 g) 

3/o(l + e)(t-U)2(i+-) 


(84) 


(85) 


In Eig. ([3]) we plot the EoS parameter as a function of time, for e = ^, fo = 0.0001(sec)“^'^“^ and tg = 10“^®sec. As 
we can see, the EoS describes a phantom evolution, for t < tg, at t tg it develops a singularity at some time during 
the phantom era, and for t > tg slowly by slowly goes from the phantom era to the de Sitter evolution. Therefore for 
t ^ U, and t <C U, the evolution is nearly de Sitter. This singular behavior in the deep phantom era has also been 
observed in other bouncing models too, see for example and we shall briefly discuss this issue later on. We can 



t(sec) 


FIG. 3: The equation of state corresponding to the cosmological bounce a{t) = ^ ^ \ as a function of time, for e = 

/o = 0 . 0001 ( 300 )“^*^“^ and ts = 10“®®sec. 


also study the behavior of the EoS analytically in various limits of the cosmic time. Let us start with the case that 
the bouncing point, and simultaneously the point that the Type IV singularity occurs, is very small, that is tg/t <C 1. 
In such a case, for large t values, and specifically for all cosmic times that t ^ tg, the EoS would become, 

= 3/0(1+ £)(t)2(l+e)’ (86) 

which for large t describes de Sitter acceleration that slightly crosses the phantom divide. Moreover, for times near 
the Type IV singularity, that is f ~ ts, the EoS becomes singular, as it can be seen from Eq. (1851) . This can also be 
seen in Eig. |31 Now consider the case that tg ^ t, in which case, when t tg, the same singular behavior occurs and 
the EoS evolves rapidly in the deep phantom era. For cosmic times much smaller than tg, that is tg ^ t, the EoS is 
approximately equal to. 


WeS = -1 - 


(1 + 2 £) 


3/0(1+ £)ts^'+"^’ 


(87) 


which again describes nearly de Sitter but slightly crosses to the phantom state. It is worth examining another 
interesting scenario, in which the EoS evolves to the singular phantom state after our present epoch, which is roughly 
tp ~ lO^^sec. In that case the EoS up to our era is nearly de Sitter, and evolves to phantom near the Type IV 


^ Equation IMli is valid when the F[R) geometric contribution is viewed as a perfect fluid, see Appendix A for details 
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FIG. 4: The equation of state corresponding to the cosmological bounce a{t) = ^ \ as a function of time, for e = ^, 

/o = 10“®^(sec)“^®“^ (left), /o = 10“^® (right) and ts = lO^^sec. 


TABLE I: Exact Valnes of the EoS for /o = 10 ®^(sec) ^ 

t (sec) _ Weft 

t = 5 X 10^® Weff = -1.000000000000058 
t = 6.1 X 10^® Weff = -1.000000000000099 
t = 7 X 10^® Weff = -1.000000000000175 
t = 8 X 10^® Weff = -1.000000000000425 
t = 1.1 X 10^^ Weff = -1.000000000001927 
f = 1.5 X 10^® Wefi = -1.000000000000058 


singularity. For example let tg = lO^^sec. Then, for t < tg, the EoS is of nearly de Sitter type and as t —>■ fg, the EoS 
becomes rapidly phantom. This can also be seen in Eig. |4j left panel, where we used /o = 10“^®(sec)“^^“^. In the 
right panel we have plotted the EoS for /q = 10“^®(sec)“^®“^. Observe that the “throat” of the EoS graph becomes 
smaller, as the value of /o increases. So practically, according to the left plot scenario of Eig. 21 this scenario suggests 
that the EoS can possibly and simultaneously cross the phantom divide, with a functional behavior that resembles to 
a Dirac function. This seems somehow like a simultaneous and peculiarly sudden change in the state of the Universe 
as a whole. It would be interesting to model such a behavior with a Dirac delta function at the transition, to see the 
differences between the bounce model ([ 3 ]) and the Dirac delta function model, but we defer this task to a future work. 
In the plots of Fig. 21 there appear many 1”, and this is due to the fact that the values of the EoS are very close 
to —1, with the differences appearing after 12 decimal places. In order to have a more clear picture of this, in TableU 
we present the exact values of the EoS for various times and for fo = 10“^^ (sec)where the differences between 
various points can explicitly be seen. Eurthermore, in Table HIl we present the corresponding exact EoS values for 
fo = 10“^®(sec)“^®“^. In order to have a clear picture of the behavior of the EoS, in Fig. Owe plot the behavior of the 
function (weff + 1) x 10^"^ as a function of the cosmic time for fo = 10“^^(sec)“^'^“^ (left) and fo = 10 “^®(sec)“^^“^ 
(right). In these plots, the differences of the EoS in the two cases can be clearly seen. The vertical line in the plots 
of Fig. O corresponds to the value t = 10^®sec. 

Let us now in brief compare the behavior of the EoS corresponding to the singular bounce (21) to the one corre¬ 
sponding to the bounce (2]). The two models as we now explicitly demonstrate are qualitatively similar, with the only 


TABLE II: Exact Values of the EoS for fo = 10"®® (sec)"®®"® 
t (sec) _ WeB 

t = 5x 10®® Weff = -1.000000000057529 
t = 6.1 X 10®® Weff = -1.000000000012679 
t = 7x 10®® Weff = -1.000000000175356 
f = 8 X 10®® Weff = -1.000000000424738 
t = 1.1 X 10®® Weff = -1.000000001927141 
f = 1.2 X 10®® WeB = -1.000000000424738 
t = 1.5 X 10®® WeB = -1.000000000057529 
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FIG. 5; The function (Weft + 1) x 10^^ corresponding to the cosmological bounce a{t) = ^ , as a function of time, 

for £ = /o = 10“®^(sec)“^^“^ (left), /o = 10“^®(sec)“^®“^ (right) and ta = 10^®sec. 


difference being that the two models are generated by different F{R) gravities in the Jordan frame. For an account 
on the bounce (l5|), see M- For the bounce ([5]), the EoS reads, 


Weff = -1 - 


1 


( 88 ) 


In Fig. [n]we plotted the time dependence of the bounce (O, for /o 


O.OOOl(sec) ^ and ts = 0. As is obvious from 



I(sec) 


FIG. 6: The equation of state corresponding to the cosmological bounce a{t) = , as a function of time, for /o = 

0.0001(sec)“^^“^ and ta = Osec. 

Fig. m and as we already mentioned, the qualitative behavior of the bounces ([3]) and ®, are very much alike. It is 
of particular importance to notice the appearance of a singularity in the EoS of the two bounces and also that both 
bounces in the limit asymptotically approach a de Sitter evolution. Moreover, near the bouncing point the two 

bounces generate a phantom EoS. Notice however that the bounce of Eq. ([5]) is not singular at all (the derivative of 
the Hubble rate is constant and all higher derivatives are simply zero). But what does a singularity in the equation of 
state indicates? In the next section we shall briefly discuss this issue and by using some illustrative examples, always 
in the context of bouncing cosmologies, we try to shed some light on this issue. 


A. Brief discussion on singularities in the effective equation of state in general bouncing cosmologies 


In this section we shall discuss the issue of singularities of the EoS in the context of bouncing cosmology. Apart 
from the singular bounce which we studied in this paper, there exist also other bouncing cosmologies, such as the 
superbounce and the matter bounce [T^ scenarios. It is worth discussing briefly how the EoS behaves in 

these bouncing cosmologies. We start off with the matter bounce scenario [l2|, for which the scale factor reads, 


a{t) = 



and the corresponding Hubble rate is equal to, 


H{t) 


2tpc 

2 + Pc 


(89) 


(90) 
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where pc is the critical energy density. This parameter frequently occurs in the LQC matter bounce scenario, and 
it can take various values depending on the model under study (see for example page 17 of [Slj, for a detailed 
discussion). The matter bounce cosmology described by the Eqs. (IMll and (IMll has all the bouncing cosmology 
features we described earlier, as can also be easily verified by looking Fig. where we plotted the time dependence 
of the scale factor a{t) and of the Hubble rate Hit), as a function of time, for pc = 10®J/m^. Of course, in the matter 




FIG. 7: The scale factor a{t) (left plot) and the Hubble rate (right plot) as a function of the cosmic time t, for ts = 10 ®®sec, 
and pc = 10®J/m® for a{t) = (fpct^ + l)® 


bounce case, the bouncing point is t = 0. The EoS corresponding to matter bounce scenario reads. 


2H{t) _ 2 

3i7(t)2 “ ~W^c' 


(91) 


Clearly, the EoS (ED is singular at the bouncing point t = 0. The case of the matter bounce scenario is very similar 
to the singular bounce case we studied in this paper, with scale factor given in Eq. m and also to the bounce of 
Eq. (E)) studied in [l^. In addition, the quintom bounce studied in Ref. @ has exactly the same characteristics as 
the singular bounce we studied in this article and also as the matter bounce scenario has, although in Q, standard 
Einstein Hilbert gravity was employed, with regards to the F{R) content of the theory. Notice that in all three 
aforementioned cases, the EoS is singular at the bouncing point, which seems to be a repeating pattern, common in 
certain types of bouncing cosmologies. We need to stress that a singular EoS in astrophysical systems is a rather 
peculiar feature which never occurs, so in principle, the same would be expected in cosmological theories. So what 
we need to understand is what does that singularity in the EoS indicates and how it is related to the finite time 
singularities we presented in a previous section. Also we have to note that in both the matter bounce and singular 
bounce case, the Hubble rate vanishes at the bouncing point, so probably a singularity in the EoS is unavoidable, 
unless this is somehow cancelled if behaves in the same way as H. However, for the functional forms of the Hubble 
rate in matter bounce and the singular bounce cases, this requirement is not fulfilled. In addition, note that the fact 
that EoS parameter rceff is singular does not always imply that the effective pressure and/or effective energy-density 
are singular too [1^, . Both of them maybe regular while the effective pressure and effective energy density might 

be finite. 

Eor example, there exist some other bouncing cosmologies in the literature, such as the superbounce [l[lS[Il, for 
which the requirement that behaves in the same way as H occurs. In this case, no singularity occurs in the EoS. 
It is worth recalling this case in brief for completeness. In the superbounce scenario [J, [^, the scale factor is 

equal to. 


and the corresponding Hubble rate is equal to, 


a{t) = {t- tg)'^ , 

(92) 


(93) 


where c is a dimensionless parameter, which is assumed to satisfy c > in order for the superbounce to occur(see [T^ 
and references therein). The superbounce (|M1) clearly describes a bouncing cosmology, having most of the properties 
of a bounce we described in a previous section, with the only exception being that at the bouncing point t = tg, 
the Hubble rate is not zero but diverges. This can also be verified by looking Eig. dZD, where we plotted the time 
dependence of the scale factor and of the Hubble rate, for c = and tg = 10“^®sec. As we can see in Fig. ®, 
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FIG. 8: The scale factor a{t) (left plot) and the Hubble rate (right plot) as a function of the cosmic time t, for ts = 10 ®®sec, 
and c = %/? for a{t) = {t — tg)'^ 


the scale factor is never singular and the Universe contracts for t < ts and expands for t > ts, where at t = tg no 
singularity occurs. In addition, the EoS for the superbounce (IMl) is equal to, 

WeS = (-3 + c^), (94) 

which is constant. So in the superbounce case, only the Hubble rate is singular, with the EoS being non-singular 
and constant and this occurs because and H behave in the same way as functions of the cosmic time. This is 
in contrast to the singular bounce we studied in this paper (jS]), the bounce (©, which was studied in [l^ and also 
to the quintom bounce which was studied in Q. Notice that in all the three aforementioned cases, at the bouncing 
point t = ts, the Hubble rate vanishes, that is. Hits) = 0. Someone could claim that for example, the singularity 
in the EoS, is due to the reason that Hits) = 0, which is true in both cases of the bounces. So provisionally, we 
could claim that the singularity in the EoS is, at a first glance, not directly related to the finite time cosmological 
singularities in the case that the Type IV singularity is taken into account and only that. This however should be 
studied in more detail and caution and therefore, the need for the complete understanding of the relation between 
the finite time singularities and the EoS singularities, is compelling. This however exceeds the purposes of this paper 
and we hope to address this issue concretely in a future work. 


VII. CONCLUSIONS 

In this paper we investigated which FiR) gravity can generate a Type IV singular bounce, which was chosen so that 
the Type IV singularity occurs at the bouncing point. As we explicitly demonstrated, the F(R) gravity responsible 
for the Type IV bounce, near the singularity has the form F{R) = R + aR? -|- A in the Jordan frame. We also found 
the Einstein frame scalar theory counterpart corresponding to the F{R) gravity and also, having assumed a quasi de 
Sitter metric in the Einstein frame, we investigated the inflationary properties of the resulting scalar theory in the 
Einstein frame. As we evinced, the spectral indices of the resulting scalar theory can be compatible with the recent 
Planck [53 observations and the theory it self is compatible with the standard Einstein frame R? inflation theory, 
with the only difference being that the minimum of the potential is shifted. 

We also performed a thorough analysis of the EoS corresponding to the Type IV singular bounce. As we showed, 
the behavior of the singular bounce is similar to other bounces that exist in the literature, such as the matter bounce 
scenario [l^ or the quintom bounce scenario Q. At this point it is worth discussing a quite interesting point that 
resulted from our analysis. As we demonstrated, even in the case of the Type IV singular bounce, the EoS is singular 
at the bouncing point, where the singularity occurs, with the other two aforementioned bounces also having this 
feature. In addition, as can be easily checked, in the case of a Type I singularity, the EoS can in some cases be 
singular or regular at the point where the singularity occurs and the same applies for the Type H case. Furthermore, 
it has pointed out in the literature that there exist points at which the EoS diverges but both the effective energy 
density and the effective pressure are finite so this would, in some sense, have nothing to do with finite time 
singularities. As for the Type H case, the EoS is always singular. It would therefore be interesting to address the 
question, whether the finite time singularities should be further classified by also taking into account the EoS. This is 
also motivated by the fact that in most known perfect fluids that are used in astrophysics, the EoS is never singular, 
so the appearance of this singular behavior in some cases, in a cosmological context should be noticed and used in 
some classifying way. This is also supported by the fact that, even in the Big Rip case, which is a crushing type 
singularity, the EoS can be regular, even though all other quantities diverge. In any case, the singularity in the EoS 
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of perfect fluids used in cosmology should be further studied in order to understand the implications this would have 
at a theoretical and at an observational level. For a relevant study on these issues, see [s^ . 

Another important issue that we did not address is the connection between the Type IV singularities and bouncing 
cosmology. Is it possible that a Type IV singular cosmological evolution would lead to a bouncing cosmology? And 
in addition, is the EoS always singular in the case of a Type IV singular cosmological evolution? We hope to address 
these issues formally in a future work. 
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Appendix A:EoS for F{R) modified gravity 


In this appendix, we demonstrate in detail that the effective equation of state of the F{R) modified gravity, for a 
flat FRW background, is given by Eq. (I5T)) . namely. 


Weff = -1 - 


2H{t) 


(95) 


For a detailed account on these issues, see also [s^ and references therein. We start off with the FRW equations of 
motion, which can be derived from Eq. o, and these read, 


3F\R)H^ = K^{p^ + pr) + - iHF'{R), (96) 

- 2F'{R)H = K^{pm + 4/3p^) + FF'{R) - HF'{R), 


so by assuming that no matter fluids are present, the above equations are modified as follows, 

3F'{R)H^ = {R)R- FiR)) _ 


- 2F'{R)H = FF'{R) - HF'{R), 

The equations (1^ can be rewritten as follows, 

3FI^ = pde, — 2H = Poe + Pde, 


(97) 


(98) 


where we defined the geometric energy density poE and the geometric effective pressure Pde, as follows. 


{F'{R)R-F{R)) 

POE = -^- 

Poe = FF'{R) + 2HF'{R) 


3HF'{R) + 377^ 


(99) 


In this way, the contribution of the F{R) gravity can be viewed as a perfect fluid’s one and the representation of the 
FRW equations appearing in Eq. (HB can be viewed as a mathematically equivalent perfect fluid representation of 
the ones appearing in Eq. (IMll . This is the reason why the contribution of F{R) gravity to the FRW equations is said 
to be owing to the existence of a geometric dark fluid. The geometric effective energy density poE and the geometric 
effective pressure density satisfy the continuity equation, 


Pde + 377 {poE + Pde) = 0, 


( 100 ) 


The corresponding effective equation of state for the dark fluid is defined to be equal to, 

Pde 

Wee = -, 

Pde 

which owing to Eq. dSIl) can be written as follows. 


( 101 ) 


Wee = - 


277 + 377(7)2 


( 102 ) 


which can easily be cast in the form of Eq. (IMj) . So for a flat FRW metric, the modified gravity effective equation of 
state can take the form of Eq. (IMj) , like in the ordinary Einstein-Hilbert gravity in a flat FRW background. 
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Appendix B: Analytic form of constants 

Here we present the full form of the constants appearing in Eqs. (1321) and (l33l) . The parameter A is equal to 
/-TO(2 + 2£)5(iTiTCir(2p^) /o“^Ci(2 + 2£)^ (l + 5e + 2£2)r(l + 


A = 


2(l+e) 


2 (i + =)r(-5ltfe) 


2 (l + e) 2 r ( 


1 _ l+4e 
^ 2 ( 1 +£) 


while the parameters B and C are equal to, 

/o 


H=6/o(1 + 4£)Ci( 

r 


i) 


2 (e+l) 


2 (e + l). 

(l + 2 (i^) ((2 + 2(e + 1)) - 2(1 + 2(e + 


(2(e + l))2r(^^ 


C = 


22 (£ + 1 ) + 1 

QfWjrr 


(2(£+l))3+2(e+l) 

(2(£ + 1))(2 + 2(e + l))r ( 2 ( 6 + 1 )) + 4e)^r(l 


(' 


2 ( 6 + 1 ) 


Y 1+46 A 

V 2 l^J 




4(^l + 2(e+l) + 2(e + lf 

)ri 

(1 + ifT+ry) 

1 4(l-k2(£-kl))2(l-b4e)r| 

(2+2(7+Ty) 

r| 

(iwfi) 

(1 + 2(£ + 1)) 


(103) 


(104) 
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